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Abstract 

We reveal unifying thermodynamic aspects of so different phenomena 
as the cosmological electron-positron annihilation, the evaporation of pri- 
mordial black holes with a narrow mass range, and the "deflationary" tran- 
sition from an initial de Sitter phase to a subsequent standard Friedmann- 
Lemaitre-Robertson- Walker (FLRW) behavior. 
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■ 1 Introduction 

The thermodynamics of two fluids with different temperatures represents a 
+3 • framework which is sufficiently general to apply to entirely different epochs 

of the cosmological evolution. This unifying feature may be used to establish 
surprising similarities between otherwise quite independent phenomena in the 
expanding universe. In this paper we focus on aspects of the temperature evo- 
lution during periods with decay and production of particles to demonstrate the 
universal power of the thermodynamic description. In particular, we show that 
the same simple law for the cooling rate of a fluid in the expanding universe 
governs a wide range of phenomena implying the cosmological electron-positron 
annihilation after neutrino decoupling at about 1 MeV, the evaporation of pri- 
mordial black holes (PBHs) with a narrow mass range, and the "deflationary" 
transition from an initial de Sitter stage to a subsequent FLRW period, 
equivalent to a phenomenological vacuum decay model. All these processes are 
characterized by a strong back reaction of decay and production processes on 
the thermal evolution of the universe. It is the possibility of taking into account 
this back reaction in a rather straightforward but general way, which admits an 
application to such a variety of different phenomena. 

More specifically, we shall first reproduce the factor (11/4) by which 
the temperature of the neutrino background differs from that of the photon 
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background as a consequence of electron-positron annihilation. Secondly, we 
show that the black hole temperature behavior r,„„. oc m _1 , where m,„„, is 
the black hole mass, is consistent with the general fluid temperature law for a 
PBH "fluid" , a configuration in which all its members are assumed to have the 
same mass rn (BH) . On this basis we discuss thermodynamical aspects of PBH 
evaporation. The third example is the evolution of the radiation temperature in 
a "deflationary" scenario of the early universe which implies an initial increase 
to a maximum value as a result of the production of relativistic particles out of 
a decaying vacuum, followed by a decrease which finally approaches the familiar 
FLRW behavior. 

None of these results is really new. The first case is cosmological textbook 
physics (see, e.g. ||), the second one was investigated in ||, the scenario char- 
acterizing the third case is based on f§ (see also |, |). What is new however, 
is the unifying view which allows the discussion of so different cosmological 
effects starting from the same set of basic equations. To highlight the under- 
lying common thermodynamical features of the mentioned phenomena is the 
main purpose of this paper. The material is organized as follows. In section 
[2] we recall the basic relations of two-fluid thermodynamics in an expanding 
universe. These relations are used in section || to discuss the cosmological 
electron-positron annihilation. In section |] the general formalism is applied to 
a mixture of radiation and a component of PBHs which are assumed to have 
the same mass. It may be shown that under this condition they share essential 
properties with a pressureless gas. Thermodynamic aspects of a smooth transi- 
tion from a de Sitter stage to a radiation dominated FLRW phase including an 
intermediate temperature maximum are investigated on the same general basis 
in section [5], while the final section ^ is devoted to concluding remarks. Units 
have been chosen so that c = £;£ = fo = l. 



2 Basic thermodynamic relations 

We assume the cosmic medium to consist of two components which share the 
same 4- velocity u l . Each of the components has a perfect fluid structure with 

the energy- momentum tensor T lk , where A = 1,2, and a corresponding particle 

{A) 

flow vector N*., 

(A) ' 

T S) = P W ^ k +P W h ik , Nl A) = n (A) u l , (A = 1,2) . (1) 

Here, p {A) is the energy density of component A, measured by a comoving ob- 
server, p {A) is the corresponding equilibrium pressure, h lk = g lk + u l u k is the 
spatial projection tensor, and n (A) is the number density of species- A parti- 



"nfe nor jyi 

(A) (A, 

interparticle conversions are admitted: 



cles. Neither T,. s nor TV* are required to be conserved, i.e., interactions and 

(A) (A) 



T l% k = ~ f U ■ = V) + 3 ^V) = n iA^ (A) ■ (2) 

The quantity H is the Hubble parameter H = a/a with the scale factor a of 
the Robertson- Walker metric. T (A) = N (A) /N (A) denotes the rate of change of 
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the number N, = n,.,a 3 of particles in a comoving volume a 3 . The T ifc and 

(A) (A) o (^4) 

AT* add up to the corresponding quantities for the medium as a whole: 

(1) ^ (2) ' JV JV (l)^ JV (2) • \°) 

It is well known that in general the energy-momentum tensor T lk does not 
take the form of a perfect fluid, but will contain a non-equilibrium pressure II 

0, i i- Different from the iJJ , the overall energy-momentum tensor has to 
be conserved, which establishes a relation between i' and t\ : 

T ifc = P u i u k + (p + n) , r'J. = o f = -t* . (4) 

We do not require, however, conservation of the total particle number [TT| , 

1. e., 

]\T« = n + 3Hn = nT = n (1) T (1) + n (2) T (2) , (5) 

where n = n (1) + n (2) is the overall particle number density. Each component is 
governed by its own Gibbs equation which provides us with an expression for 
the time behaviour of the entropy per particle s. 



T (aAa) = d^+j^d-— ' V) T {V(A) = **f W -(P W +P W ) F (A) 

(A) n (A) 

(6) 

In general, the temperatures T {A) of both components are different. With the 
help of the equations of state 

P(A) =P( A) (V)' T (A)) > P W =P ( A) (V)' T (A)) ' ( ? ) 

one obtains the evolution law for the temperatures T (A) . Namely, differentiating 
p. in (Q) along the fluid flow lines and applying the balances (||) we find 

1,0,0 

7 " 3ff , (8) 



\ 3H)dp {A) dp {A) /dT {A) 



where 



d P(A) y dp ^ /dT ^;n (A) dp w fdp w 

d P(A) (dp w /BT w ) ' 8T W [dT ( 



(A) 7 n (A) 

The general temperature law @ provides the unifying basis for the discussions 
of the following sections. It will play a central role in our investigations of 
both the electron-positron annihilation and the PBH evaporation and a specific 
inflationary scenario. 

An important special case which we frequently will refer to is characterized 
by the condition s (A) = 0, which means constant entropy per particle fi"o| , |i"Tf . 
This condition simplifies the temperature law, 
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Moreover, according to (Q) it establishes the link u a t^ = (^p (A) +P (A) J 
between the source terms I\ and t l which together with the last relation of 

(A) (A) ° 

(||) provides us with a relation between the rates and T (2) : 

(2) n (2) 

where h {A) = (p {A) +P {A ^j l n ^ A ) are ^ e enthalpies per particle. 

3 Electron-positron annihilation 

Let us consider the cosmological period of electron-positron annihilation and 
the corresponding creation of photons, shortly after neutrino decoupling (see, 
e.g. |D, §3.1.2). Annihilation becomes predominant as soon as the radiation 
temperature drops below the mass of the electron. Before the electron-positron 
annihilation there are two bosonic degrees of freedom (photons), four fermionic 
degrees of freedom due to the electrons and positrons, and 12 fermionic degrees 
of freedom due to the different neutrino species. The fermionic energy density 
is 

7 vr 4 3C(3) ^3 nu 

PF= 630CW nF7> ' nF = l~^3FT F , (11) 

where £ (x) is Riemann's Zeta-function, and the bosonic one 

PB = 30 ^ ( 3 ) UbTb ' nB = ^~9bTI . (12) 

The factors gp and gs are the numbers of fermionic and bosonic degrees of 
freedom, respectively. 

After the electron-positron annihilation we are left with the two bosonic 
degrees of freedom (photons) and the 12 neutrino degrees of freedom, i.e., four 
fermionic degrees of freedom have disappeared. The neutrino degrees of freedom 
are not affected at all by this process. The neutrino temperature behaves as 
T (v) = T. , (to) a (to) /a(t), where the initial time to is assumed to be a time 
before the beginning of the annihilation process, i.e., electrons, positrons and 
photons are still at equilibrium at to with T (v) (to) = T, ± . (to) = (to) = T (0) . 
Let us consider the subsystem of four fermionic and two bosonic degrees of 



freedom [13]. The four fermionic degrees of freedom due to the electrons and 
positrons are identified with fluid 1 of our general analysis, i.e., (1) — > (e^) 
while the two bosonic degrees of freedom due to the photons are fluid 2, i.e., 
(2) — > (7). With these specifications the number densities become 

n (e±) ~ 4 4 7j-2 i (e ±) ' n (7) ~ Z n 2 1 M ■ W 



The corresponding enthalpies per particle are 

14 i 4 4 vr 4 

V) ~ YWW) {e±) ' (7) ~ 330C"(3) 
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Assuming h, ± ./h^ to be given by their (constant) initial ratio, Eq. (|lC|) 
integrates to 

iV (T) (t) = N h) (t ) + ^ [jV (e±) ( io ) - JV (a±) (t)] . (15) 

(7) 

The final value N, , (tf) corresponds to the case where the electrons and positrons 
have been annihilated, i.e., N, ± . (tf) = 0: 



N M (*/) _ , , ^(e±) V) (*0) 



1+ O ^ . (16) 

N {j) (t ) h (j) n (7) (i ) 

This result has to be coupled to the temperature law (||), which for photons 
becomes 

T (l) a 3A^ (7) w a(t) \N {i) (t )J 

For f > t f the ratio iV (i) /N, (t ) is fixed by the value (t f ) /N (t ). 
Since h {e±) /h M = 7/6 [cf. (£§] and n (e±) (t ) /n (iy) (t ) = 3/2 [cf. ©] we 

obtain iV- , (t/) /^ r (7) (*o) = 11/4. Consequently, the temperature evolution law 
for t > tf is 

T W-T °( < o)/ rll V /3 r w^_r n V /3 (t>t) Mtt 

w ( } " (o) "^y It J - It J ( * - ^ • (18) 

Thus we have reproduced the well-known difference between photon and neu- 
trino background temperatures on the basis of the temperature law (§) . Usually, 
this result is obtained by calculating the entropy transfer from the pairs to 
the photons under the condition of entropy conservation 0]. 



4 Evaporation of primordial black holes 

In a variety of scenarios with copious production of primordial black holes the 
latter may substantially contribute to the energy density of the universe (see, 
e.g., H). Some of these models are characterized by a narrow mass spectrum 
1 14]. Under such circumstances it is a good approximation to ascribe the same 
mass rn (BH) to all members of the population. On the other hand, a black 
hole mass is known to be characterized by a temperature T {BH) cx m 7 BH s • Con- 
sequently, with a single mass population of PBHs one may associate a single 
temperature T, BH . as well. Furthermore, one may show that a PBH population 
in the expanding universe may be regarded as an ensemble of non-interacting 
particles ||. These properties suggest a description of the PBH component as 
a pressureless "fluid" , in which T, BH . in some respect plays the role of a fluid 
temperature. Since PBHs are known to evaporate, it is tempting to estab- 
lish a two-fluid model along the lines of Sec. |2| with one component being the 
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PBH "fluid" , the second one radiation. The equations of state (0) for the PBH 
component are 



P 



(BH) 



^ ' P(BH) ~ n (BH) m (BH) ■ (19) 

The black hole temperature is related to its mass by the well-known formula 

1 



T, 



(BH) 



8nm 



(20) 



(BH) 



This temperature is attributed to each PBH individually, i.e., primarily it is not 
a conventional fluid temperature. The number -^ {SH) of PBHs in a comoving 
volume a 3 , N {BH) = n (BH) a 3 , is not preserved and, according to Eq. (Q), we 
may write down a balance equation for the corresponding PBH number flow 
vector N 



(BH) 



n (BH) U 



(BH), 



'"'(bh) + 3Hn {BH) 



n V 

(BH) (BH) 



(21) 



The black hole energy balance becomes [cf. Eq. (J2J) with (|l|)] 

P(BH) + ^HP(BH) = U °t\ BH) 



(BH) 



in 



(BH) 



+ 



(BH) 



in 



(BH)} 



(22) 



Using p 



(BH) 



as well as Eq. dq) in the fluid temperature law IM) we find 



u a t a , 

™ (BH) 



(BH) i (BH) 



(BH) 



in 



(BH) 



'(BH) 



d P(BH)l dT (BH) d P(BH)l dT ( 



(BH) m (BH) 



(23) 



We emphasize that we have used here the same symbol, T, gH . , for the the PBH 
"fluid" temperature and for the temperature (|20|), which is ascribed to the 
individual black holes. The consistency of this identification becomes obvious 
if we combine the equations of state (19) with ( |20| ) and introduce the result for 

d P(BH)/ dT (BH) int ° &■ 



dp, 



(BH) 



dT, 



(BH) 



(BH) 



-(BH) 



T, 



m 



(BH) 



(BH) 



'(BH) 



111 



(24) 



(BH) 



It is the crucial point of our analysis that Hawking's temperature law ( |20D for 
individual black holes fits together with the general fluid temperature law (||) 
for the equations of state { JT9| ) with ( p0| ) . This circumstance provides the basis 
for our thermodynamical discussion of the PBH evaporation process. To this 
purpose we identify component 1 of the general analysis in section ^ with the 
PBH "fluid" and component 2 with ulrarelativistic matter (radiation, subscript 
r), i.e., (1) — > (BH) and (2) — > (r). For the latter we require constant entropy 
per particle, i.e., [cf. Eq. @] 



0) 



Uat a {r) 



3^V) X (r) 



_M 
r 



■HI 



(25) 
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Combination with t a 

(BH) 



-t a from 

(r) 



yields 



r 



(r) 



4 V) 

3 



(r) 



r 



(Sir) 



+ 



(M) 



(26) 



n (BH) S (BH) U 



The total entropy flow 5 a is the sum of the contributions 

and S? = n (r) s (r) u a . With s (SH) = 47rm^ s we obtain the following expression 
for the overall entropy production density 



,1 ~~ P(BH) (BH) 



1 



1 



2T T 

^(BH) ■ L {r) 



+ P,BH 



(BH) 



m 



(BH) 



(BH) 



-(r) 



(27) 



It is obvious from (|2q ) that negative values of rh {BH) (and T (SH) ) correspond 
to a positive quantity I\ This case is equivalent to the creation of radia- 
tive particles at the expense of the PBH mass (and its number), i.e., to PBH 
evaporation. The inverse process, namely "accretion" with rn (BH) > and 



r <r) < is described by the general formula (27) as well. Which of the two 



processes is thermodynamically preferred depends on the ratio of the temper- 
atures. Given a specific initial ratio, the further evolution is entirely governed 
by the temperature laws in ( ^4|) and (25). Let's assume an initial configuration 



with T, DD ., (to) = T (to)- A non-negative entropy production density then 



(BH) 



requires r, Bfr) < 0. Since one expects Th, BH . and T {SH) to have the same sign, 
this implies a positive value of I\ , , i.e., radiation particles are produced which 
makes the PBH masses shrink. The further evolution depends on a subtle in- 
terplay between the rates T (BH) and T (r) and their respective back reactions 
on the temperature laws (24) and (]25|). A positive I\ . may either be larger 
or smaller than the expansion rate 3H. For r {r) < 3H the fluid temperature 
decreases according to Eq. (p5j), while the BH temperature increases according 
to Eq. (p4[). It follows that T {BH) > T (r) at t > to- The evaporation process 
will continue since T,_, < T, D „, requires F (BH) < and rh IBH) /m % < to 



"(r) 



guarantee 5.^, > in Eq. 



(BH 

m 



For T 



(BH) ^ (BH)/ (BH) 

> 3H, however, hypothetically real- 



ized e.g. by a large initial ratio p {BH) / p,, , the fluid temperature increases. If 
this increase is smaller than the increase in T /D „, we have again T. . < T,„„. 

(BH) ° (r) (BH) 

and the PBH evaporation goes on since it remains thermodynamically favored 
(S? a > 0). But an increase in T (r) stronger than that in T^ BH) results in a fluid 
temperature which is higher than T. BB .. For T. . > 2T (BH) a positive entropy 



production 



requires V 



(BH) 



> 



and ih, BH Jfn/ BH ) > 0, implying a quick 
transition to a negative I\ j, i.e., the process can no longer continue. A strong 
"reheating" of the fluid will stop the evaporation and reverse the process. Now, 
the second law requires PBHs to be formed out of the radiation and to accrete 
mass. A negative r (r) , on the other hand, will make T, . subsequently decrease 
[cf. Eq. (|25|)1. If T, . has fallen below T. BB ., the evaporation process may set 
in again. In particular, this self-confining property implies that a catastrophic 
growth of the PBHs is thermodynamically forbidden. The point is that a PBH 
growth, i.e. mass accretion with rii (BH) /'m {BH) > 0, back reacts on the temper- 
ature of the ambient radiation. For a fixed PBH number, i.e. r {flH) = 0, the 
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corresponding radiation temperature changes as 



T (r) = _ H ( 1 + 1 P(BH) m (BH) H -A 
T ir) _ V 4 P(r) m (B H) J ' 

It is obvious that for fh, BH , > from some time on the temperature T {r) will 
cool off faster than T (BH) [cf. Eq. (PH)]. Consequently, Tl. will approach 
T, BH y As soon as T {r) has fallen below T, BH ., mass accretion stops since for 
T (r) < T (BH) the rate m, OT ,/m (M , has to be negative in order to guarantee 
a positive entropy production, i.e., the process now proceeds in the reverse 
direction and the PBHs can no longer grow but start to evaporate again. This 
completes our thermodynamic discussion of PBH evaporation based on the 
temperature law (|8|) (and its special case @)). 



5 "Deflationary" universe 

In the two previous examples we did not consider the impact of the decay and 
production processes on the expansion behavior of the universe. As was shown 
in |l3j and [H], the general tendency of this influence is to increase the cosmic 
expansion rate. Namely, processes of the type discussed in sections [3| and |] give 
rise to an effective viscous pressure of the cosmic medium as a whole [cf. Eq. 
(H)]. Since this contribution to the overall pressure is negative, its net effect is 
to accelerate the expansion. While this effect is small for the cases dealt with 
in sections || and it is essential in the "deflationary" universe model of the 
present section. This model relies on Einstein's field equations with the energy- 
momentum tensor (Q) of a bulk viscous fluid. In a homogeneous and isotropic 
universe one has 

k P = 3H 2 , H = -^(p + p + U) kU = -3-/H 2 - 2H , (28) 

where k is Einstein's gravitational constant and 7 = 1+p/p. In case IT is not a 
"conventional" viscous pressure but represents a quantity describing cosmolog- 



ical particle production on a phenomenological level [10, pdj| , it may be related 
to the production rate T introduced in (||) . For "adiabatic" particle production 
this relation is 

n = -(p + p)^. (29) 



(30) 



Combination with the field equations fl28|) then yields || 

2 H H' 37 



1 + 



3H 3 7 H 2 H 



-E- - 1 

3H L 



2 a 



where H' = dH/da. Strictly speaking, the rate T has to be calculated on the 
quantum level (see, e.g., fi~6| , |l7|). In a phenomenological setting an ansatz for 
T/H is required. For a dependence r oc p oc H 2 jl], ||, || and 7 = 4/3 we obtain 

a 2 e T a\ 

a z + a% 3H a z + ai 



S 



where we have chosen the constants such that H e = —H^, i.e, a > for a < a e 
and a < for a > a e . H starts with a constant value Hq = 2H e at a <C a e 
and then "deflates" towards the typical H oc a~ 2 behaviour of a radiation 
dominated universe for a S> a e . This Hubble rate has originally been obtained 
in the context of phenomenological approaches to cosmological vacuum decay 
||, |J. Again, this is a two-component model with one component playing the 
role of the cosmological "vacuum" . Our point here is to demonstrate that such 
kind of model fits into the general structure of section ^ and admits a similar 
thermodynamic discussion as the cases of electron-positron annihilation and 
PBH evaporation. We will identify the first component of the general formalism 
in Sec. |2|with the "vacuum" (subscript v), i.e., (1) — > (v), the second one again 
with radiation, i.e., (2) — > (r). The sketched scenario may then be obtained on 



the basis of an interacting two-fluid model with p 



P M + P. 



(r) 



where 



>) 



3H 2 2 



i 3 



a 2 + a 2 



M 



ml 

2vr 



m 2 P 





[ *l 


(«e) 


a 2 + a 2 



1 3 



(32) 



and nip = 8tt/k is the square of the Planck mass. The part p, . is finite for 
a —>■ and decays as a~ 6 for a S> o e , while the part p (r) describes relativistic 
matter with p, > — > for a — > and p.. oc a -4 for a ^> a e . The energy balances 
are (A = v, r) 



p (A) +3H p (M +p 



'(A) 



(A) 



\A) 



+ P 



{A) 



(33) 



with 



3i? 



2 a 2 ) a 2 + a 2 



(r) 

3H 



2a 2 + a 2 



The equation for p may be written as 

p. , + 3ff ( p, , + p. . + n, , 



(34) 



(35) 



where 

n 



3H 



2 ol / a 2 + 

This corresponds to an effective equation of state 



tfiPw +p w) 



(36) 



a 2 - a 2 



(37) 



Although we have always p, . = /0 ( ,/3, the effective equation of state for a — > 
approaches P, . = ~Pr v y Effectively, this component behaves as a vacuum 
contribution. For a > a e it represents stiff matter with P.. = p... The 
radiation component may be regarded as emerging from the decay of the initial 
vacuum according to 

p (r) +AHp (r) =-p (v) . (38) 
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The radiation temperature is obtained from the general law ([]), which in the 
present case specifies to 



T. 



O) 



(r) 



-HI 



3H 



Integration yields 



T, 



w 



2 3/4 T ^£ 



a 2 + a 2 



3/4 



2 a 2 + a 2 



T, . (x p 

(r) r 



(39) 



1/1 

(r) 



(40) 



T, is the value of the radiation temperature at a = a e . This temperature 



' ( e . r ) 

starts at T, 

r 



(r) 



— for a = 0, then increases to a maximum value, given by 
3H, equivalent to a 2 - 1 



rpmax 

(r) 



32X 1 / 4 
27/ 



(e,r) 



(41) 



and finally decreases as a -1 for large values of a. Our formalism allows us 
to ascribe a temperature T,-. to the "vacuum" as well, which is determined 
analogously by 



r 

"(») 



-HI 



3fT 



The "vacuum" temperature behaves as 



-(«) 



a 2 + a 2 



2^a 2 + a 2 ' 



3/4 



(42) 



(43) 



It starts from a maximum value at a = and decreases as a 3 / 2 for large a. The 

"vacuum" is radiative in the sense that p. , oc T 4 is valid. As a final remark 

r M (u) 

we mention that it is also possible to introduce a temperature T of the cosmic 
medium whole with a behavior || 



r 

3^ 



-H 



a 2 + o 2 



.1/2 



a 2 + a 2 



(44) 



which "interpolates" between (p3D for small a and (pC|) for a S> a e - These 
considerations clarify the central role played by the general temperature law 
(^) also under conditions where the relevant back reaction substantially affects 
the entire cosmological dynamics. 



6 Conclusions 

Cosmological thermodynamics allows us to establish a unifying view on a broad 
range of different phenomena and to uncover joint underlying structures. In this 
paper we have explored similar thermodynamic features of matter creation in 
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the early universe, primordial black hole evaporation, and electron-positron 
annihilation after neutrino decoupling. All these processes are governed by the 
same basic temperature law for a fluid with variable particle number which takes 
into account the back reaction of the relevant interactions on the thermal history 
of the universe. A particular aspect of our considerations is the consistency of 
this law with Hawking's black hole temperature formula T. D „, oc . This 

circumstance provides the basis for a two-fluid model for the evaporation of a 
single-mass PBH component into radiation. 
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